Macroscopic quantum coherence in spinor condensates confined in an anisotropic 

potential 
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We investigate the macroscopic quantum coherence of a spin-1 Rb condensate confined in an 
anisotropic potential. Under the single-mode approximation, we show that the system can be 
modeled as a biaxial quantum magnet due to the interplay between the magnetic dipole-dipole 
interaction and the anisotropic potential. By applying a magnetic field along the hard-axis, we 
show that the tunneling splitting oscillates as a function of the field strength. We also propose an 
experimental scheme to detect the oscillatory behavior of the tunneling splitting by employing the 
Landau-Zener tunneling. 
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I. INTRODUCTION 

Tunneling of a macroscopic variable into a classically 
forbidden region provides one of the most striking mani- 
festations of quantum mechanics [1| . The model of quan- 
tum tunneling often involves a particle moving in a mul- 
tistable potential. Quantum mechanically, there is a fi- 
nite probability for the particle to tunnel through the 
barrier and escape from a metastable state to an ab- 
solutely stable one, which is often referred to as the 
macroscopic quantum tunneling (MQT). In a symmet- 
ric double-well, the particle tunnels through the barrier 
to oscillate back and forth between the degenerate states, 
known as macroscopic quantum coherence (MQC). The 
quantum tunneling removes the degeneracy of the ground 
states, resulting a tunneling splitting between the true 
ground state and the first excited state. 

Because of its small size and precise characterizability, 
single-molecule nanomagnet represents an ideal platform 
for demonstrating the MQT and MQC of the spin 
For such a system, MQT consists of the tunneling of the 
magnetization out of the metastable easy directions in 
the presence of external field, which was experimentally 
observed as a series of steps in the hysteresis loops in 
Mni2-acetate and Fes molecules at low temperatures @- 
9]. MQC in spin systems represents the resonance be- 
tween two equivalent easy directions. Of particular in- 
terest, in the present of a general biaxial anisotropy, the 
spin has two preferred tunneling path via the medium 
axis when the external field is applied along the hard 
direction. As a result, the constructive and destructive 
interference then give rise to the oscillatory tunneling 
splitting [10(. The geometric-phase effect of the quan- 
tum spin tunneling was first observed by Wernsdorfer 
and Sessoli (Tlj . 

In the context of the ultracold atomic gases, the MQC 
problem was previously studied for a spinor conden- 
sate trapped in a double-well potential, where the mag- 



netic dipole-dipole interaction (MDDI) between the con- 
densates confined in different wells induces an uniaxial 
anisotropy with the easy-axis b eing along the direction 
connecting two potential wells [12j. In fact, even in a 
single axially symmetric trap, a dipolar spinor conden- 
sate can be treated as an uniaxial quantum magnet [l3j |. 
whose magnetic properties, spin squeezing, magnetiza- 
tion steps, and macroscopic entanglement generation, 
and MQC were studied [3 EH- However, for uniaxial 
quantum magnet, the tunneling splitting of model is a 
monotonically increasing function of the transverse field 
strength [l2j ■ The possibility of obtaining the oscillatory 
tunneling splitting in ultracold atomic gases was explored 
by considering a condensate coupled dispersively with an 
ultrahigh-finesse optical cavity [la ]. 



In the present work, we investigate the MQC of a 
spin-1 Rb condensate confined in a three-dimensional 
anisotropic harmonic oscillator potential. Under the 
SMA, we showed that the interplay of the MDDI and 
the anistropic trap results in a biaxial quantum magnet 
whose magnetic structure can be approximately specified 
by the geometry of the trapping potential. Subsequently, 
we study the MQC of of the condensate by applying an 
external magnetic field along the hard-axis. Similar to 
the molecular magnet, we show that the tunneling split- 
ting of our system oscillates as a function of the field 
strength. Finally, we propose an experimental scheme 
to detect the MQC by utilizing the Landau-Zener tran- 
sition. 



This paper is organized as follows. In Sec. [Til under the 
single-mode approximation (SMA) , we derive the Hamil- 
tonian of the system. In Sec. lIIIl we explore the magnetic 
structure of the system. Section IIVI is devoted to the 
properties of the MQC and its experimental detection. 
Finally, we conclude in Sec. [V] 
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II. MODEL act via the s-wave collisions and the MDDI. In the sec- 

ond quantized form, the total Hamiltonian of the system 
We consider a trapped gas of TV spin F = 1 Rb atoms reads [l3| 
subjected to an external magnetic field B. Atoms inter- 
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where ip a is the field operator for mF — a spin state, 
V^ xt is the anisotropic confining potential, F is the an- 
gular momentum operator, gp is the Lande g- factor, 
and fiB is the Bohr magneton. The spin-independent 
and spin-exchange s-wave collsions are characterized by 
c = 47rn 2 (a + 2a 2 )/(3Af) and c 2 = 47rft 2 (a 2 -a )/(3M), 
respectively, with a/=o,2 being the scattering length of 
two spin-1 atoms in the combined symmetric channel of 
the total spin / [T3, [H|]. In particular, we have c 2 < 
for Rb atoms, indicating that the spin-exchange inter- 
action is ferromagnetic. The strength of the MDDI is 
Cd = MoA*b.9f / '(47r) with fj,Q being the vacuum permeabil- 
ity. 

To proceed further, we adopt the SMA which assumes 
that atoms in different spin states share a common spatial 
mode function <ft(r). The field operators can then be 
decomposed into [l9( 



ip a (r) ~ <j){r)a c 



(2) 



where a a being the annihilation operator of the spin com- 
ponent a. The validity of the SMA can be justified by 
noting that Co 3> |c 2 | and Cd — 0.1|c 2 | for Rb atom |20| . 
By substituting Eq. into ([1]) and dropping the con- 
stant terms, the total Hamiltonian reduces to 

H = (c 2 - c^S 2 + 3c' d S 2 z - 34(S 2 X - S 2 V ) - g Fl i B B ■ S 



+ 3c' d ald Q + 3c^'(aLiai + af 



(3) 



where S = J2 a p ^Fapap is the total many-body an- 
gular momentum operator and S v (r) = x, y, z) is its 
projection along 77-axis, c 2 = (c 2 /2) J dr|</>(r)| 4 is the 
strength of spin-exchange interaction, and the strength 
of the MDDI arc characterized by two parameters: c' d — 
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r'l J sin 2 1? e 2 ^ 



with "d and ip being the polar and azimuthal angles of 
the vector r — r', respectively. It can be shown that 
c'j = if the mode function <fi(r) possesses an axial sym- 
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FIG. 1: (Color online). Anisotropic constants D(K x ,K y ) (a) 
and E(k x , k v ) (b) for Rb condensates. 



([3]) originates from the commutation relations between 
the bosonic operators. 

It is convenient to rescale the Eq. (j3)) by using \c' 2 \ as 
the energy unit, which yields the dimensionless Hamilto- 
nian 
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metry [13j . One should note that the last line of the Eq. 



where D = —3c' d /\c' 2 \ and E = — 3c^'/|c 2 | are the 
anisotropy constants and H = gpfisB /\c' 2 \ is the dimen- 
sionless magnetic field. Apparently, the Hamiltonian ([5]) 
describes a biaxial quantum magnet. To be more spe- 
cific, we assume that the mode function is a Gaussian 
0(r) = 7r- 3 /% x g yg2 )-V2 e - WK) with % be _ 

ing the width of the condensate along the rj direction. It 
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FIG. 2: Magnetic structure of a Rb condensate on the pa- 
rameter (k x , K y ) plane in terms of (easy, medium, hard) axes 
in the absence of magnetic field. The solid lines correspond 
to uniaxial magnets. 



In the absence of the external magnetic field, the mag- 
netic properties of the reduced Hamiltonian ([7]) are com- 
pletely determined by the anisotropic constants D and 
E. In Fig. [2 we classify the magnetic structure of the 
system on the geometric parameter plane in terms of the 
easy, medium, and hard axes, which essentially states 
that the easy (hard) axis corresponds to the direction 
with the weakest (strongest) confinement. Figure can 
be physically understood as follows. Without the MDDI, 
the ground state energy is degenerate about the orienta- 
tion of the spin. The MDDI removes this degeneracy by 
pointing the spins to certain directions. Noting that, for 
two magnetic dipoles, the dipolar interaction is attractive 
(repulsive) for the head-to-tail (side-by-side) configura- 
tion. Therefore, the spins prefer to align in the way such 
that the possibility for the head-to-tail configuration is 
maximized, which is exactly the direction corresponding 
to the weakest confinement. 
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47rcd 
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where the set of parameters (k x ,k v ) = (q x /q Z j Qy/Qz) 
characterizes the geometry of the condensate, io.i(') is 
the modified Bessel functions of the first kind, and erfc(-) 
is the complementary error function. Figure Q] shows the 
values of D(k x , k v ) and E(k x , k v ) for Rb condensate. In 
particular, E{k x , k v ) — when k x = 
becomes an uniaxial magnet [l3|, E3] • 



K y , the condensate 



III. MAGNETIC STRUCTURE 

To find the ground state wave function, we have nu- 
merically diagonalized the Hamiltonian Q in the Fock 
state basis, {\Ni,No,N-i)} subjected to the constrains 
N a > and J2 a = N, for the parameter regime cov- 
ering N < 50, 0.1 < K x<y < 10, and arbitrary magnetic 
field. It is found that, as long as N > 5, we always have 
(S 2 ) ~ N(N + 1), indicating that the condensate stays 
in a ferromagnetic state with total spin S ~ N. In ad- 
dition, the contribution from T-L' [Eq. ©] is negligible. 
Consequently, we may drop the constant S 2 term in Hq 
and the linear term T-L' such that the Eq. (Q| reduces to 
the familiar biaxial Hamiltonian 

H cS ~ -DSl + E{S 2 X — Sy) — H • S (7) 

for quantum magnet with total spin S — N . 



IV. MACROSCOPIC QUANTUM COHERENCE 

A simple model for studying the MQC problem is an 
uniaxial magnet with easy-axis anisotropy, which can be 
realized, for example, by taking D < and E = in the 
Hamiltonian For ultracold atomic gases, the MQC of 
such model was previously studied in the Refs. [l2l. fl5| . 
Here, without loss of generality, we study the MQC of a 
biaxial magnet with D > E > 0, which corresponds to 
the geometric parameters (k x ,k v ) in the shaded region 
of Fig. [2j The magnetic field is applied along the hard 
axis, i.e., H = H x ~k. 

We first consider the classical counterpart of the re- 
duced Hamiltonian (J7]) by treating S as a vector of length 
|S| — N. The total energy then becomes 

£(6, 4>)=- DN 2 cos 2 6 + EN 2 sin 2 6 cos 2(f> 
— H X N s'm9 cos <fi, 

where (6, </>) denotes the direction of S. Apparently, 
for H x < H* = 2N(D + E), there exist two degener- 
ate ground states located at (#o,0) and (n — 9 ,0) with 
sin6» = H X /[2N(D + E)]. When H x > H*, the dou- 
ble degeneracy is removed such that the system is fully 
polarized along the x axis by the external field. 

Quantum mechanically, the classical degeneracy of the 
ground states is lifted by quantum tunneling even when 
H x < H*, which results a tunneling splitting A£ be- 
tween the true ground state and the first excited state. 
Of particular interest, with the biaxial anisotropy, there 
exist two tunneling paths with opposite windings on the 
yz easy anisotropy plane. The constructive and destruc- 
tive interferences of quantum spin phases of the different 
paths cause the tunneling splitting to oscillate with the 
magnetic field. This phenomenon was first predicted by 
Garg [Io| and was experimentally observed by Werns- 
dorfer and Sessoli in Mni2 molecules [11]. Using the in- 



stanton method, the tunneling splitting can be expressed 




analytically as [10|, |2l| 

AS = Ae |cos(7r6)|, 



(8) 



where @{H X ) = N - H x / 2^/2E(D + E) is the area on 

the Bloch sphere enclosed by the two instanton paths and 
Aeo is the tunneling splitting under zero external field 
which contains the contributions from the classical action 
and the fluctuations around the instanton paths [T3, [H| ■ 
Clearly, the quantum tunneling is completely quenched 
whenever = n + 1/2 with n being an integer. The 
period of this oscillation is [l(| 



AH X = 2y/2E(D + E). 



(9) 



We want to emphasize that the instanton method is only 
valid for H x < H* as it depends on the classical paths of 
the tunneling. 

In Fig. [31 we present the field dependence of the tun- 
neling splitting obtained via the exact numerical diag- 
onalization of the full Hamiltonian ((4]) for N = 10, 
lg k x — —0.9, and \gn y = —0.1. As a comparison, we 
also plot Eq. © with Ae from Ref. [2l| . As can be seen, 
for small H Xl the agreement between two approaches is 
quite well, which further confirms that our system are 
well described by the reduced Hamiltonian. However, 
significant discrepancy is observed for the larger external 
field. It is found that, by increasing the atom number 
N, the agreement between the numerical and instanton 
methods can be improved. 

We now turn to discuss the experimental detection of 
the MQC in spin-1 Rb condensates. Here we adopt a sim- 
ilar scheme to that used in the molecular magnet experi- 
ment by utilizing the Landau-Zener transition ll|. More 



specifically, the experiment can be carried out as follows. 
In addition to the constant transverse field iJ x x, we intro- 
duce a time-dependent longitudinal magnetic field H z (t)z 
which is swept linearly with a constant rate v > 0, i.e., 



FIG. 4: (Color online). No/N as a function of H x for tp = 
(a) and tp with H z ^ — 0.2 (b). Other parameters are 
N = 10, lgn x = -0.9, \gK, y = -0.1, v = 10 -3 , and = 



H z (t) = H [ z ] + vt for t > 0. When the longitudinal field 

reaches designated = + vtf > at t = tf, one 
measures the number of atoms N a in each spin compo- 
nent. 

As can be seen from the reduced Hamiltonian (J7]), if the 

initial longitudinal field (< 0) is sufficiently large, 
the initial state roughly stays at the S z = —N level which 
anticrosses with the S z = N level at H z = 0. When the 
longitudinal field sweeps through this avoid crossing, a 
Landau-Zener transition occurs. Since the energy gap 
between the levels S z = —N and N at the anticrossing is 
exactly the tunneling splitting AS which depends on the 
transverse field, it is expected that N a will oscillate as a 
function of H x . 

The proposed experiment can be simulated by numer- 
ically evolving the full Hamiltonian ((4]) with an initial 
state being the ground state under the initial magnetic 
field. Figure H)Ja) shows the typical transverse field de- 
pendence of No for various final longitudinal fields. As 
expected, Ao oscillates as a function of the transverse 
field strength. The oscillation period is also in very good 
agreement with that of AS. In addition, by increasing 

(f) 

H z , more anticrossings in the energy spectrum of the 
Hamiltonian are swept through. As a result, the oscil- 
lation amplitude decreases with , and eventually the 
Nq(H x ) will roughly converge to the curve corresponding 

( f) (f) 
to m j ' = 8 as one further increases H z . Experimen- 
tally, it is possible that the transverse field is misaligned 
such that it forms an angle ip to the a;-axis. In Fig.0Jb), 
we plot Nq(H z ) for tp 0. The oscillation remains for 
small if, however, it will disappear for large tp. Finally, 
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we remark that N± and 7V_i also exhibit the similar os- 
cillatory behavior as N (H Z ). 

V. CONCLUSION 

To conclude, we have studied the MQC of a spin-1 
Rb condensate confined in an anisotropic trap. Under 
the SMA, we showed that this system can be described 
as a biaxial quantum magnet. Physically, this biaxial 
anisotropy is induced by the MDDI and the anisotropic 
trap. Subsequently, we showed that the magnetic struc- 
ture of the system is determined by the geometric pa- 
rameters of the trapping potential. We then studied the 
MQC of the spinor condensate by applying an external 



magnetic along the hard-axis, we showed that the tunnel- 
ing splitting oscillates as a function of the field strength. 
Finally, we propose to detect the MQC by utilizing the 
Landau-Zcncr transition. 
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